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Abstract 

We study the exclusive production of tttt and ptr in hard 7*7 scattering in the 
forward kinematical region where the virtuality of one photon provides us with a 
hard scale in the process. The newly introduced concept of Transition Distribution 
Amplitudes (TDA) is used to perform a QCD calculation of these reactions thanks 
to two simple models for TDAs. Cross sections for pit and tttt production are eval- 
uated and compared to the possible background from the Bremsstrahlung process. 
This picture may be tested at intense electron-positron colliders such as CLEO and 
B factories. The cross section e7 — > e'ii 7i is finally shown to provide a possible 
determination of the ir° axial form factor, F\ , at small t, which seems not to be 
measurable elsewhere. 



1 Introduction 

In a recent paper [T] , we have advocated that factorisation theorems j2] for exclusive processes 
may be extended to the case of the reaction 7r~ 7r + — > 7* 7 in the kinematical regime where 
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the virtual photon is highly virtual (of the order of the energy squared of the reaction) but 
the momentum transfer t is small. This enlarges the successful description of deep exclusive 
reactions in terms of distribution amplitudes j3j and/or generalised parton distributions 01 E] 
on the one side and perturbatively calculable coefficient functions describing hard scattering 
at the partonic level on the other side. We want here to describe along the same lines the 
(crossed) reactions 7^ 7 — > AB: 

7z,7 ~ ^ P ±7r=F ; itl^ if^ir*, 7^7^vr 7r°, (1) 

in the near forward region and for large virtual photon invariant mass Q, which may be 
studied in detail at intense electron colliders such as those which are mostly used as B 
factories. 



7*(g) ^ A{p A ) 




Figure 1: The factorised amplitude for 7*7 — > An at small transfer momentum. 

Let us recall the main ingredients of the analysis developed in Ref. QjJ Ej. With the 
kinematics described in Fig. Q and introducing light-cone coordinates v ± = (v° ± v 3 ) / \f2 
and transverse components vt = (f 1 ,^ 2 ) for any four-vector v, we define the 7 — > n transi- 
tion distribution amplitudes (TDAs) T(x,£,t) as the Fourier transform of matrix elements 
(vrQ^IO^G^)) where CP = $(-§)[-§ , f] r^(| ) with V 1 = t^V,^" We then 
factorise the amplitude of the process 7^7 — > An as 3 

M{Q\ (,t) oc J dxdz $ A (z)M h (z, x, C)T(x, (, t) , (2) 

with a hard amplitude Mh(z,x,£) and &a(z) is the hadron distribution amplitude (DA). 

The variable z is as usual the light-cone momentum fraction carried by the quark entering 
the meson A, x+l; (resp. x— £) is the corresponding one for the quark leaving (resp. entering) 
the TDA. The skewness variable £ describes the loss of light-cone momentum of the incident 
photon and is connected to the Bjorken variable xb (see section 3 for detailed kinematics). 

Contrarily to the case of generalised parton distributions (GPD) where the forward limit 
is related to the conventional parton distributions measured in the deep inelastic scattering 

3 The 7^ case is more difficult to analyse since the leading twist amplitude vanishes. 
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(DIS), there is no such interesting constraints for the new TDAs. The only constraints are 
sum rules obtained by taking the local limit of the corresponding operators and possibly 
soft limits when the produced-meson momentum vanishes. Lacking any non-perturbative 
calculations of matrix element defining TDAs we are forced to build toy models to get 
estimates for the cross sections, to be compared with future experimental data. In particular, 
Tiburzi [7 has recently developed a model for both these axial and vector TDAs, which we 
shall use in a simplified version to make a comparison with the results obtained from our 
model. Also some other approaches used to model GPD in the pion case |H] could be applied 
in this context. 



2 The 



7 



7T 



TDAs 



Let us first stress an obvious point. The 7 



7T J 



TDA involve of course only quark correla- 



tors, and so does the 7 



7T 



TDA since the charge conjugation property of the ir° forbids 



the leading twist gluonic TDA. Therefore, in the following, we only need to take into account 
quark correlators. For definiteness, let us consider the 7 — > tt~ TDAs which are given by PQ 
(P = S ^ L , A=p ir --j h ): 
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(3) 



where the first two TDAs, V n (x,£,t) and A n (x,£,t) (commonly labelled by G(x,£,t) in 
the following) are chiral even and the latter ones, T[ (x,£,t),i = 1, 2, are chiral odd. In the 
TDA Eq. (JHJ), we include the Wilson line [y; z]=P exp [ig(y - z) J^dt n^ity + (1 - t)z)}, 
which provides the QCD-gauge invariance for non local operators and equals unity in a 
light-like (axial) gauge. On the other hand, we do not write the electromagnetic Wilson line 
caused by the presence of the photon, since we choose an electromagnetic axial gauge for 
the photon. 

The TDAs V n (x, £, t), A* (x, £, t), (x, t) and T% (x, t) are dimensionless quanti- 
ties; f n is the pion decay constant (/„. = 131 MeV). These four leading twist TDAs are in fact 
linear combinations of the four independent helicity amplitudes for the process 57 — ► qn~ . 
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2.1 Constraints on 7 — > 7r TDAs 



Sum rules may be derived for the photon to meson TDAs. Since the local matrix elements 
appear in radiative weak decays, we can relate them to the vector and axial form factors Fy 
and Fa- 

For 7r ± , for which one has [§| 

= 0.017 ± 0.008 and Ff = 0.0116 ± 0.0016, (4) 

we get 4 

/ dxG**(x,Z,t) = —Ft(t) G=(V,A). (5) 

To what concerns the neutral 7r°, one can constrain the vector TDA thanks to the elec- 
tromagnetic (transition) form factor F 7r o 7 * 7 which is known at small t j^j- Indeed, one has 5 

J dx (Q u Vf (x, £, t) + Q d Vf (x, £, tj) = U F^ 7 (t), (6) 

where Q u = 2/3, Q d = —1/3 and Vg°(x,£,t) is the TDA related to the operator built 
from the quark q. Current algebra fixes the value of the right hand side at t — since 

^-r*-r(* = 0) = DH- 

These sum rules 6 constrain possible parametrisations of the TDAs. Note, in particular, 
the ^-independence of the right hand side of the relations. 

2.2 Toy models for 7 — > it axial and vector TDA 

Since we want to respect the £ independence of the first moment in x of the TDA leading to 
the sum rules (Eq. (JHJ) and Eq. @), we shall use the analogy with the construction of GPDs 
through the double distributions [TT]. We consider here the case of the 7 — > n~ TDA. 

2.2.1 t-independent double distributions 

In this section, we start from double distribution for t — [11 and we omit the D-term which 
being isoscalar does not exist for the 7 —>■ TDA. A D-term may exist for the 7 — >• 7r° 
TDA, but we shall not include it for simplicity. 

We start from the representation of the x and £ dependence of the TDA in the form 

G^(x,0= f dp f m d a 6(x-P-Sa)f(/3,a), (7) 
J-i J-i+\P\ 

4 In Eq. JSJ), we use the same definition as in PDG 9 , Fq being a dimensionless quantity. According to 
these conventions, the Fq defined in should be divided by m n . 

5 This sum rule is slightly different from the one written in pQ where we had assumed a quark model 
decomposition of the tt°. 

6 In principle, one can also derive analogous sum rules for the chiral-odd TDAs (T\ and T2 in Eq. |J3J|) but 
we do not see an obvious way to relate them to experimentally interesting observables, see the discussion at 
the end of section [3] 
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with 



f(P,a) = q(P)h(P,a) 



where q(/3) is analogous to the forward quark distribution in the GPD case. h(/3, a) is a 
profile function usually parametrised as 



h^(P, 



a) 



T(26 + 2) [(1 



2 -a 2 ] b 



2 26+i r2 ( 6+1 ) (l-l/^l) 



2b+l 



(9) 



where the parameter b characterises the strength of the £— dependence. As a first guess we 
assume that the (3— dependence of q is given by a simple linear law 



q{(5) =2(1-/3) 9(0) . 



(10) 



As usually done for GPDs, we assume a mild £ dependence as given by b = 1 and implement 
the normalisation (with J dxG^(x,^) = 1) and t— dependence of the TDA through the axial 
or vector form factor: 



G(x,U) = G^\x,0-—F G (t) 



The t— dependence of these form factors has been studied in chiral perturbation theory ^3] 
and turned out to be weak, so we shall neglect it in this model and we shall use the measured 
values of Eq. (J1J) . 
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+ (1 - 0(f - x))2(x - 1) 2£ + lo. 



1 + ^ 



(12) 



We do not include any QCD evolution to our TDAs, the effects of which are supposedly 
much less important than the uncertainty of our modelling. The resulting 7 — > 7r~ vector 
TDA is plotted on Fig. |21 (a). The axial one and the ones for 7r° have the same shape with 
a different normalisation. In the following, we shall refer to this approach as Model 1. 
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2.2.2 t-dependent double distributions 

We shall use here the initial t-dependent double distribution of Tiburzi Explicitly, we 
have: 



with 



v(x,z,t) = ^ I dp '' 



ni 



da S(x - (3 - £a)W(P, a; t), 



i+ 



A(x,£,t) = — d(3 da6(x-P-Za)B(/3,a]t), 

m 7T J -I J-1+\P\ 



(13) 
(14) 



W(P,a;t) 
B(P,a;t) 



m 



47T 

3m 2 {a + P) 



t 

'4 



m 2 + -P(a + 0- \)m\ - [I - a 2 - P{2 - P)] 

m 2 + -P{a + P- l)m 2 - [1 - a 2 - P{2 - p)\ 



An 2 



(15) 



(16) 



where m is set to 0.18 GeV [7|. The sum rule for V is satisfied for t = —0.5 GeV 2 . We shall 
keep t fixed in the following to enable comparison with Model 1. In the following, we shall 
refer to this approach as Model 2. We see that these two models give quite different results 
and we shall use both in section 3 to estimate the sensitivity of the cross sections to TDA 
models. 



V(x,U) 



V(x,^,t) 





Figure 2: The 7 — ► tt vector transition distribution amplitude V(x,£,t) in Model 1 and in 
Model 2 (for t = -0.5 GeV 2 ). 
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3 e7 — > ep^7r cross-section estimates 



We first consider the pi vr production case and more precisely when the p flies in the direction 
of the virtual photon and the 7r enters the TDA. It is cleaner than the 7r ± 7r T production 
case. Indeed, the reaction 7*7 — > AB is accessible only in a e7 reaction 7 and thus must 
compete with the Bremsstrahlung process where the meson pair An comes from a virtual 
photon emitted by the lepton line. In the p^ir case, this Bremsstrahlung process is much 
suppressed because of the vector and axial- vector character of the pi, and n leading-twist 
DAs; this leads to a vanishing leading-twist contribution to the form factor F p7T at high 
transfer momentum. Moreover, in the neutral p° 7r° case, the TDA process is forbidden by 
C-conjugation. 



3.1 Kinematics 

The momenta in the process ej — > ep ± n :r are defined as shown on Fig. El in the CMS of the 
p and 7r T (or equally of the 7 and 7*). We choose the x — z plane to be the one of the 
collisions 7*7 — > pit, which we call hadronic plane. The leptonic one, i.e. where the 7* is 
emitted by the electron is at an angle ip of the hadronic one. We also have P = P7rT 2 +P7 , 
A = p n ^ — p 7 and q = p e — p' e , q being the momentum of the 7*. 




Figure 3: The kinematics of e(p e ) + 7(p 7 ) — > e(p' e ) + p(p p ) + n(Pw) i n the center of mass of 
the meson pair. 

Then we define the light-cone vectors p and n (p 2 = n 2 = 0) such that 2p.n = 1 and the 
invariants in terms of scalar products of the momenta: 

Q 2 = -q 2 = -{Pe-p' e )\ £ = 

* = (PttT - Py) 2 = A 2 , S n = (p e + p y ) 2 , 

^ - (n _L t, ^ „ _ 1-P-i _ (g+l)Q 2 _ Q 2 + W 2 



W 2 = (q + Pj ) 2 , y 



P-y-Pe 2£s e 



7 We do not consider the complete case of e + e collisions which may as usual be rewritten in the equivalent 
photon approximation in terms of the e 7 reaction and the quasi real photon flux. 
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For definiteriess, we choose, in the CMS of the meson pair, p = ^^^ (1, 0, 0, — 1) and 
n = 2 (Q2+w 2 ) ^' ^' ^) an< ^ we can then express the momenta trough their Sudakov decom- 



position 8 : 



2_Lj4/2 ,n2 



Pr = (1 + OP, q = ^Stn - + £)p 



p?r- = (i-0p-^ + a t , At = ^ 

We can see that £ is determined by the external kinematics of 7*7 — > p7r through £ ~ 
q2^2vf 2 ~~ similarly to = q-^w 2 ^° wm ch it is linked via the simple relation £ ~ 2 ^ . 

Since we want to focus on the study of the TDA behaviour, we decide to choose Q 2 , t, £ 
and tp as our kinematical variables. The differential cross section thus reads 

J e-y-^epir i i 

■ITW^^I 2 . (17) 



dQ 2 dtd£d<p 32(2vr) 4 s2 f (f + 1 



3.2 The 7^7 — > 7r T amplitude 

Let us first consider the longitudinally polarised p meson, described by its twist-2 DA. Now 
the vector character of this DA selects the vector TDA and the amplitude for the reaction 
7l7 ~^ Pl n ~ becomes proportional to V(x, £, t). It reads 



M™ A (Q 2 ,£,t) = - I dx f dz^<f> p (z)M h (z,x,£)V(x,£,t) 



'-1 JO 

where the hard amplitude is (z — 1 — z) 

M h{ z, z , i) = i^p^± ( + «t_A ^ n ^ vAc . (19) 

N C Q zz\x — £,+ie x + £ — ze / 

In the hard subprocess, the quark momenta have been as usual restricted to their com- 
ponent collinear w.r.t the associated meson, i.e. for the quark leaving (resp. entering) the 
TDA, k ~ {x + £)p (resp. k' ~ {x — £)p) and the quark (resp. antiquark) entering the p 
meson, £ ~ zp p (resp. £' ~ zp p ). 

We choose <j) p (z) = 6zz as the asymptotic normalised meson distribution amplitude, 
f p = 0.216 GeV [Hj. After separating the real and imaginary part of the amplitude, the 
^-integration gives: 

4 V = t dx( Q ; t , + Q * . )v(x,t,t) 



1 \x — £ -Me x + £ — it 

Qu f 1 ix ^,(,t)-nu,t) +Qi /■ dg v(».{.o-v(-t,t,«) 

7-1 y_! x + £ 

+g^(e,e,t)(io g -i7r)+Q d y(-e,e,*)(iog f^|J +^) 



(20) 



8 A| < 0. 
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The scaling law for the amplitude is 

M™ A (Q 2 ,t,t) 



Q 



(21) 



up to logarithmic corrections due to the anomalous dimension of the TDA. 
3.3 The 7^7 — ► p% tt t contribution to cj —> ep L ir T 

The squared amplitude for cy — > epf n T from the subprocess 7^7 — > pf n T is obtained 
from J\A^I^ A (Q 2 , £,) and from the contribution of the fermionic line for the emission of a 
longitudinal virtual photon 



\M 



TDA 1 2 

e7 



Ana 



2Q 



(22) 
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Figure 4: ^-dependence of the ratio of the real part to the modulus of the TDA amplitude 
for t = -0.5 GeV 2 and Q 2 = 4 GeV 2 . 

Averaging over the real photon polarisation and integrating over f thanks to the p- 
independence of the TDA process, we eventually obtain the differential cross section: 



da 



TDA 

e7 — >ep+7r~ 



1287raLa 2 27T fl 



-t){2s n i - (£ + 1)Q 2 )(1 - £)(Re 2 (/J) + Im 2 (^)). (23) 



dQ 2 dtdi 9(£ + l) 4 Q 8 s e7 /2 



First, we see that the sole dependence on t is through the factor —t which comes from 
e^^n^EuPpArj and physically accounts for the penalty to turn a helicity ±1 state into 
a helicity 0; indeed the cross section vanishes at t = 0. Secondly the positivity of the 
kinematical factor (2s e7 £ — (£ + l)Q 2 ) yields an upper bound for Q 2 for given s n and £. The 
s ei dependence of the cross section comes entirely from the factor ( 2se ^~(£ +1 )Q ) _ 

Both real and imaginary part of the amplitude contribute significantly to the cross section, 
which is reasonable at these moderate energies. To quantify this statement, we plot on Fig.|3] 
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the relative contribution of the real part of the amplitude, as a function of the skewness £. 
It is independent of the other variables. 

The phenomenological analysis of the pion form factor ^5] seems to indicate that a rather 
large (~ 1) value of a s should be used together with the asymptotic DA. In our plots, we 
therefore put a s — 1. Our upcoming conclusions would not be strongly affected by a different 
choice. 




Figure 5: cj — > e'p^ir differential cross section plotted as a function of £ (a) and Q 2 (b) for 
s ei = 40 GeV 2 , t = -0.5 GeV 2 and respectively Q 2 = 4 GeV 2 for (a) and f = 0.2 for (b). 

In Fig. |3] (a), we plot the cross section vs £ and in Fig. |3] (b) vs Q 2 for both models of the 
vector TDA. The ^-dependence can be generically understood at large £ by the factor (1 — £) 
and at small values by the limitation on the phase space by the kinematical factor previously 
discussed. The behaviour for intermediate values of £ is sensitive to specific models for the 
TDA. As shown in Fig. 0(b), the Q 2 -behaviour is model independent and thus constitutes 
a crucial test of the validity of our approach. 

Let us finally add a remark about the case of transversally polarised p, which could have 
been quite interesting since one may naively try to write the amplitude of the TDA under 
the form 

M p *;(Q 2 ,0oc J dxdz<S>f> T (z)M h (z,x,OT(x,Z,t) , (24) 

with the chiral-odd DA for jEj and the tensor chiral-odd TDAs since chiral-odd quantities 
should appear in pairs in a physical amplitude. However, a straightforward calculation shows 
that such an amplitude vanishes at leading order due to the identity r y a iT fJiUr y a = 0. This is 
reminiscent of the analysis of transversally-polarised vector-meson electroproduction in the 
forward region, where it has been shown ,fBJ that chiral-odd GPD contributions vanish at 
all orders and the measurement of these GPDs needs another hard subprocess |17j . 
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4 e7 —> en n cross-section estimates 

4.1 7^7 —> 7T + 7r amplitude 

The momenta in the process cy — > eTr + ir~ are defined as in the c-f —> ep + tt~ case discussed 
above, where the 7r + meson flies in the direction of the virtual photon and the 7r~ meson in 
the direction of the real one. 

The main difference with the previous case is that the pion DA replaces the p one and that 
the axial 7 — > 11 TDA A(x,£,t) defined in Eq. (jHJ) replaces the vector one. The amplitude 
thus reads : 



1 rl 

1 

1 



M™ A {Q 2 ,£)= I dx I dz(j) n (z)M h (z,x,OA(x,Z,t), (25) 



where the hard amplitude is 

/ - 4:n 2 a em a s C F 1 / Q u Q d \ 

M h {z,x,£) = — — r 7—— H 7-7 Je-A, 26 

N c Q z z \x — 4 + it x + 4 — 16 J 

and where <fi(z) is the normalised meson distribution amplitude, z = 1 — z. The factor 
/tt cancels with the one from the TDA definition. Now if we choose <p°^{z) = 6zz, the 
^-integration is readily carried out and the result has the form of Eq. (|2U|) with obvious 
replacement of TDAs V — > A. The scaling law for the amplitude is similar to the one for 
the p7r case. 



4.2 The e"f —> e7r + tt cross section 

As said above, the contribution of the subprocess 7^7 — > 7r + 7r~ to the ej — > eii + tt~ cross 
section is obtained from A1™ A (Q 2 ,£) with the addition of the fermionic line for the emis- 
sion of a longitudinal virtual photon. The cross section is obtained from Eq. (j23jl by the 
substitutions I v — > J A and / p — > 

The Bremsstrahlung process where the 7r + 7r~ pair is produced by a photon radiated from 
the leptonic line (see the graphs of Fig. EJ) has the following squared amplitude (omitting for 
now the interference terms) 



M 



cy — >e , n~ r K 



with 



.4 
B 
C 
D 



|2 64(4vr« em ) 3 e 2 



-t 



-t t 



m + \^\^+^ + ^B + ,l~C + ^D) (27) 



Q 2 Q 



Q 1 



COS(f- 



2(3£-l)(y-2) 



[£(5£ - 2) + l]y 2 - 2(y - !)[(! - 3£) 2 + 4(£ - l)£cos2^] 



COS<£>- 



4C(3e-l)(y-2) 
1 + 0vT+IvT^?vT^1/ 



K + i) 
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Figure 6: Feynman graphs for the Bremsstrahlung contribution to — > e7r + 7r at LO. 

The factor used at the 7r + 7r~ vertex is ie(p n - — p 7T +)' 1 F n ((q + p 7 ) 2 ) where F n (W 2 ) is the 
pion form factor which is well measured for W 2 up to a few GeV 2 and is described in a 
satisfactorily way at large space-like values by perturbative QCD [3 , while the time-like 
region is less understood [THj but is measured [2T)] . 

In the following estimates, and to be consistent with our previous choice of the asymptotic 
pion DA in the TDA subprocess amplitude, we describe the pion form factor F^iW 2 ) with 
the asymptotic distribution amplitude and a large value of a s . This choice will not affect 
the main conclusions of this section, and the inclusion of a fit to the available data would 
give the same order of magnitude for the Bremsstrahlung contribution, which will remain 
mostly negligible in all the kinematical domain considered. 

The use of the asymptotic form for the pion DA has been criticised and other inputs 
have been proposed on the basis of QCD sum rules [21]. If one uses the C-Z parametrisation 
and choose a s = 0.6 , the amplitude remains unchanged. This has to be paralleled with the 
form factor analysis, where the use of the C-Z DA is sometimes accompanied by the use of 
a lower value for the strong coupling constant a s of the order of 0.4 |15j . 

We show in Fig.[7|the relative contributions of the Bremsstrahlung and TDA subprocesses 
(omitting interference terms) to the differential cross section integrated over ip, as a function 
of s e7 , £ and Q 2 . Except for small values of s ei , the TDA subprocess clearly dominates for all 
reasonable values of £. We also show the 92-dependence of the Bremsstrahlung contribution, 
which may be turned into a positive check of the dominance of the <p independent TDA 
subprocess, much in the manner of the successful test of the GPD framework in deeply 
virtual Compton scattering 2'2\. 

Since the Bremsstrahlung contribution is quite small, we have not given the expression 
of the interference terms which may be anyhow cancelled out if one considers the charge- 
averaged cross section da('K + iT~) + <ier(7r~7r + ) since the TDA and Bremsstrahlung amplitudes 
produce hadronic states with opposite charge conjugations. On the other hand, this latter 
property may be used to separate the interference contribution by measuring the charge 
asymmetric quantity da(n + 7i^) — <io~(7r~7r + ) , and thus analysing the TDA subprocess at the 
amplitude level [23J. 
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2 3 4 5 6 7 -1 -0.5 0.5 1 



Q 2 coscp 

Figure 7: Differential cross sections for the TDA and the Bremsstrahlung subprocesses as a 
function of s e7 , £ and Q 2 , and for the sole Bremsstrahlung subprocess as a function of cosip 
(lower right plot). Except for the variable studied, Q 2 = 4 GeV 2 , £ = 0.2, t = —0.5 GeV 2 
and s e ^ = 40 GeV 2 . 



4.3 The neutral case: e^j — » en it 

In this case, C invariance forbids any contribution from the Bremsstrahlung subprocess. We 
emphasise here that the measurement of the cross section for cy — > eir° tt° in the forward 
region would provide one of the sole possible determination of the axial form factor of 
the 7T° at small t. Indeed, its measurement is not feasible in 7r° — > vv^j and it is completely 
drown in the electromagnetic background in 7r° — > e + e~7. 

Therefore, any measurement of the cross section for e7 — > e7r° 7r° would constrain the 
value of F\ via the sum rule for the axial TDAs A^°(x, £, t) similar to the one linking F„-o 7 * 7 
to Vg°(x, £, t) (see Eq. (jHJ)). We expect the cross section 7*7 — > 7r° 7r° to be of the order of 
the one for 7*7 — > n T depicted in Fig. 
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5 Conclusion 



In summary, we have shown how to test experimentally the new factorised QCD approach 
to forward hard exclusive scattering. We believe that our models for the photon to meson 
transition distribution amplitudes are sufficiently constrained to give reasonable orders of 
magnitude for the estimated cross sections. The good news is that the hard hadronic process 
dominates the Bremsstrahlung contribution in the kinematics which are accessible in existing 
e + e~ colliders. The goal is now to test our approach, in particular by verifying the scaling 
of the cross sections, and then to measure these new hadronic matrix elements. Cross 
sections are large enough for quantitative studies to be performed. On the other hand, 
lattice studies could calculate these TDAs, at least within some approximations. Chiral 
perturbation theory may even be used to improve the extrapolation from a large quark 
mass to a realistic model. Higher-order corrections to the hard scattering process should be 
studied. The long history of the improvements of the QCD understanding of form factors at 
large momentum transfers [24] should give us a way to better include the effects of e.g. end 
point regions and their partial Sudakov suppressions. As in all studies of this type, possible 
higher- twist contributions may be relevant at measurable values of Q 2 . Some studies have 
been done in related cases in particular models j2Sj- Experimental verification of the scaling 
laws are the first tests to be applied to our description. Moreover, the kinematical domain 
of applicability of the present description (and in particular the small transfer, large energy, 
requirements) may not be clear cut in real experiments using existing accelerators. An 
interesting further study should be to try to understand the transition regions between the 
kinematics where our framework should apply on the one hand and on the other hand the 
other existing QCD descriptions of the same 7*7 — > AB reactions at small energy |2Hj or at 
very large energy [2*7] . 
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